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ABSTRACT 

With the advent of high-resolution infrared spectrographs, Radial Velocity (RV) searches enter into a new domain. As of today, the 
most important technical question to address is which wavelength reference is the most suitable for high-precision RV measurements. 
In this work we explore the usage of atmospheric absorption features. We make use of CRIRES data on two programs and three 
different targets. We re-analyze the data of the TW Hya campaign, reaching a dispersion of about 6 m/s on the RV standard in a time 
scale of roughly 1 week. We confirm the presence of a low-amplitude RV signal on TW Hya itself, roughly 3 times smaller than the 
one reported at visible wavelengths. We present RV measurements of Gl 86 as well, showing that our approach is capable of detecting 
the signal induced by a planet and correctly quantifying it. 

Our data show that CRIRES is capable of reaching a RV precision of less than 10 m/s in a time-scale of one week. The limitations of 
this particular approach are discussed, and the limiting factors on RV precision in the IR in a general way. The implications of this 
work on the design of future dedicated IR spectrographs are addressed as well. 

Key words. Planetary systems, Infrared: stars, Instrumentation: spectrographs, Methods: observational, Techniques: radial velocities 



1. Introduction 



The discovery of a Hot Jupiter orbiting 5 1 Peg by Mayor & 
Queloz (1995) triggered the quest for extrasolar planets. This 



campaign made use of the radial velocity technique (henceforth 
RV), by which the presence of a planetary companion is inferred 
by the wobble it induces on the parent star. Several dedicated sur- 
veys have employed this method, concentrating in the beginning 
on solar-type stars, the most favorable target for high-precision 
measurements. However, soon it became clear that extrasolar 
planets were ubiquitous around these hosts. Longing for a better 
understanding of formation processes, planetary searches started 
to diversify themselves, targeting both young (Soz zetti et al.| 
| 2006||Setiawan et al.|2007| > and evolved ([Lovis & Mayor|2007l 
|Sato et al.|2008) , l ow ma ss (Bonfi l7eTai.|2004||Endl et al [2008) 
and massive stars (jHatzes et al. 2005; Galla ndet al.|2005j L 
Recent technological developments allowed for more pre- 



cise spectrographs to be built, such as HARPS ( |Mayor et al. 
|2003) , and reduction and analysis methods perfected themselves 



throughout the years ( |Lovis & Pepe|2007| >. As of today, HARPS 
yields the most accurate RV measurements, with sub-m/s preci- 
sion, allowing for a succession of ground-breaking detections of 
the lightest planets known ( |Lovis et al.|2 006 ; Ma yor et al.|2 009 ). 

Both the search for more exquisite stellar hosts and the race 
for the lightest planets point towards M dwarfs. These stars, the 
most abundant in the Universe, are also the lightest ones. Since 
the RV variation induced by a planet on a star scales with M~ 2 ^ 3 , 
the amplitude of the effect induced on an M star is significantly 



* Based on observations taken at the VLT (Paranal), under programs 
280.C-5064(A) and 60.A-9051(A), and with the CORALIE spectro- 
graph at the Euler Swiss telescope (La Silla). 



larger. As an example, a planet of identical mass at the same 
distance from the stellar host produces an RV variation with an 
amplitude ~ 3 times larger on an M5 star than on a G2 star. The 
drawback is that M dwarfs, being much colder, are much fainter 
in optical wavelengths. The RV survey of light-mass stars points 
then towards the exploration of a new wavelength domain, the 
IR, where the luminosity of these objects peaks. 

As RV searches started to target young and active stars, an- 
other advantage of the IR became apparent. RV signals can be 
created by surface inhomogeneities, like stellar spots ( |Queloz| 
et al. 200 1 ] > . Since the nature of these phenomena can bypass 



current diagnosis method such as bisector analysis (see, for in- 
stance, Desort et al. 2007|>, a planetary origin can be assigned to a 



periodic variation which is, in fact, of stellar origin. By observ- 
ing in the IR, as in opposition to the visible wavelengths, one 
observes in a domain where the contrast between stellar spots 
and the stellar disk is reduced. Since the spot/ star contrast ratio 
is at the root of the RV signal detected, the amplitude of the lat- 
ter is reduced as well ( |Martm et al.|20 06 ; Huela mo et alJ 2008 ; 
|Prato et al.|2008[ ). By comparing RV drawn from optical and IR 
spectrographs, one can assign the planetary or stellar nature to a 
RV signal in an unambiguous way. 

However, the exploration of the RV in the IR is in its infancy, 
and many problems remain to be tackled. In this work we con- 
centrate ourselves in describing a precise wavelength calibration 
technique for the IR domain. In section 2 we discuss the difficul- 
ties of wavelength calibration in the IR and present the principles 
of our technique. In section 3 we describe CRIRES and the ob- 
servations, and how we applied our approach. Section 4 presents 
an outline of the data reduction and analysis. Section 5 presents 
our results. Section 6 discusses these results and their implica- 
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tions on the design of spectrographs dedicated to RV measure- 
ments in the IR. We finish by stating our conclusions in Section 
7. 



2. Wavelength References in the IR 

In order to detect small-amplitude RV signals (i.e., of m/s) a si- 
multaneous wavelength reference must be employed. In the vis- 
ible wavelengths, the Th-Ar emission lamps and the I2 absorp- 
tion cell have proved to yield a precision of better than 1 m/s 
( [Mayor et aT1|2009| ) and 3 m/s ( |Howard et~aT1|2009| ), respec- 
tively. Impressive as they are, the transfer of this knowledge to 
the IR domain is not straightforward. The density of Th-Ar lines 
is much lower in the IR than in the visible, making its usage dif- 
ficult. On top of that, since the lamp lines do not follow the same 
instrumental profile (IP) variations than the science targets (they 
pass through different optical paths), a stable spectrograph is re- 
quired. The more frequent gas-cell approach does not provide 
for a better solution, as no known gas creates a dense absorption 
forest in a wide wavelength range in the IR. Several hybrid alter- 
natives were put forward to attempt to solve this problem, such 
as gas cell plus emission lamp or coupling of different gas cells. 
For a detailed description of the problem and proposed solutions 
the reader is referred to the work of Mahadev an & Ge| ( |2009| l. 

A different approach can be undertaken by using one of the 
peculiarities of the IR: the absorption of light by the Earth's 
atmosphere. Our atmosphere is composed mainly of molecular 
species, whose rotational and vibrational features correspond to 
near-IR wavelengths. Usually seen as a nuisance, they can be 
used to our advantage. In some narrow wavelength domains, the 
presence of well-defined, sharp, and deep lines can provide for 
a good local wavelength solution. While this is only possible lo- 
cally, it is well suited for the current high-resolution IR spectro- 
graphs, characterized by a small simultaneous wavelength cov- 
erage, typically of /1/50-70. As a free-of-charge, always present 
gas cell, the usage of the atmospheric features as a wavelength 
reference is to be assessed before venturing into other alterna- 
tives. 

However, our "gas-cell" suffers from a particular problem: 
the lines it imprints on the spectra depend on the atmospheric 
conditions. While some of lines' properties are expected to re- 
main constant, others, such as depth, vary at a several % level. 
It is then very hard to apply the typical deconvolution technique 
( Butler et al.|19 96), which assumes a detailed knowledge of the 
absorbing spectra to calculate the IP and correct for it. While 
several dedicated models have been developed ( [Bailey et al.| 
2007) to address this problem, there is no evaluation of the im- 
pact of the remaining residuals on RV, and is arguably bigger 
than the m/s level. As a matter of fact, the only campaign that 
employed the deconvolution approach in the IR in the more stan- 
dard way, Blake et al. (2007 1, reported a precision of 300 m/s 
and stated that a precision down to 100 m/s was within reach. 
However, this result is to be taken with caution, as it was derived 
for a campaign of M dwarfs, an exigent stellar host. 

We decided to depart from this approach, by avoiding the 
overlapping between target and reference features. While the 
spectra of the star and of the atmosphere are simultaneous in 
the sense that they are obtained at the same moment, the lines are 
well separated spectroscopically and can be fitted independently. 
The center of each line is thus determined in a straightforward 
way (for instance, the atmospheric lines' depth variation has no 
impact on the measured RVs). The IP variations that affect stellar 
lines affect reference lines as well. If the same fitting procedures 
are employed on both line sets, the effect cancels out. 



Several notable experiments on the usage of atmospheric ab- 
sorption features were made in the eighties, of which we high- 
light the results of |Balthasar eFaX] (19821 ), |Smith| ( [T982| , and 
|Caccin et al.| ( 1985] >. Using O2 lines, these authors could reach a 
precision of 5 m/s. The value is made even more relevant by the 
fact that they used different RV determination methods (different 
instrumentation, different line fitting approaches, etc.), now be- 
lieved to be obsolete when compared with the Cross Correlation 
Function (Bar anne et al. 1996[ ) or the Deconvolution method. 
Very recently Snellen (2004) used the same principle on H2O 
lines and reached a precision of 5-10 m/s on UVES data. This 
shows that the accuracy of our wavelength reference is at least 
of the order of 5 m/s, and a valuable option for us. 

3. Observations 

3.1. CRIRES 

CRIRES (CRyogenic InfraRed Echelle Spectrograph) is the 
new ESO IR spectrograph mounted on the Unit Telescope 1 
(UT1, Antu) of the European Southern Observatory Very Large 
Telescope ( Kaeufl et al.|2004| ). The main optical elements con- 
sist of a prism acting as a pre-disperser and a 31.6 lines/mm 
echelle grating. The instrument provides for a resolution of up to 
100 000 when used with a 0.2" slit. It operates in the range 960- 
5200 nm with an instantaneous wavelength coverage of ~ A/ 50. 
The spectra is imaged on a detector mosaic, composed of four 
Aladdin III detectors (4096 x 512 pixel) with a gap of ~250 pix 
between the chips. As a consequence, the instantaneous wave- 
length coverage is much smaller than for optical cross-dispersed 
spectrographs, such as HARPS, HIRES or UCLES. Adaptive 
Optics (MACAO - Multi-Applications Curvature Adaptive op- 
tics) can be used to optimize the signal-to-noise ratio and the 
spatial resolution. 



3.2. The observing set-up 

At the time of observations, CRIRES had 200 possible observing 
set-ups, most of which overlapped. Q In order to select the most 
suitable set-up for precise RV we need to decide which one car- 
ries the higher content of useful spectral information. We stress 
here that by "useful" we mean domains where both reference - 
i.e. atmospheric - and target - i.e. stellar - spectral features are 
not superposed. 

We used solar NSO/Kitt Peak FTS data ( |Livingston &] 



Wallace 



199T] |Wallace et a"L][T993) to simulate CRIRES spec- 



tra. To do so we convolved the spectra with an IP which we as- 
sumed as gaussian and defined as a function of resolution (fixed 
at R = 100000), and re-sampled it on the detector. This proce- 
dure was applied for the Sun spectra, atmospheric spectra, and 
the product of the two. In order to access the RV uncertainty 
measured on a sun-like spectrum using the atmosphere as refer- 
ence we calculated the final RV uncertainty as the squared sum 
of the estimated error on the atmospheric lines' RV and on the 
stellar lines' RV. It is important to stress that the objective was 
not to obtain an estimation of the value of attainable precision 
in absolute terms but in a relative way between the different set- 
tings. 

To estimate the precision on each of the components two dif- 
ferent methods were employed: 1) a simple rule of thumb that 



1 Since May 2009, there are 282 standard wavelength setups avail- 
able. The standard wavelength setups are designed to provide a com- 
plete coverage with no contamination by adjacent orders. In addition, 
users can select a free wavelength setting. 
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states that the inverse of RV accuracy e oc S/N . /. yfh~i, where 
S/N is the Signal-to-noise ratio of the spectra, I is the average 
depth of absorption lines and n; the number of lines in the spec- 
tra; 2) the Quality factor described in Bouchy et al. (2001 ) which 
quantifies the spectral information present on each spectra. 

It is clear that the obtained precision depends heavily on the 
useful spectral domain. In our particular case, the superposition 
of telluric and stellar spectra can deem several zones as use- 
less for RV calculation. In order to take into account this effect, 
we proceeded as following. We used an automated algorithm to 
identify stellar lines in the simulated spectra ( Sousa et al.|2007| l. 
From these we selected only stellar lines which were not super- 
imposed on atmospheric lines bigger than a fixed threshold, and 
vice-versa. We added a Doppler shift in the interval [-30,+30] 
km/s in order to reproduce the line shift created by Earth's move- 
ment around the Sun (for the case of a star which has an average 
RV of relative to us). Lines that were found to be superposed 
after this step were discarded as well. We varied the detection 
and superposition thresholds at the level of 1-5%, a reasonable 
value for high S/N spectra. As the superposition of a relatively 
small line on the measured one can already introduce a RV shift 
of several hundreds of meters per second, it was mandatory to 
keep the superposition threshold low. 

At the end of this procedure two zones were identified as the 
most promising: 

- the H band, at 1 .6 ymi - The telluric absorption is dominated 
by CO2 lines, while other absorbers have negligible absorp- 
tion, smaller than 1%. The CO2 spectrum is easily identi- 
fied and provides for very regular wavelength anchors, with 
around 30 deep lines - from 10 to 60 %- on several detectors. 
The stellar spectra provides deep lines, which, while not be- 
ing very numerous, allow for an average contrast of around 
30%. 

- the K band, at 2.32 jjm - The atmospheric absorption shows 
features of several molecules, the most prominent being of 
CH4. The stellar lines depth is below 10%, but the high num- 
ber (around 50) compensates for it. This was the setting cho- 
sen by |Blake et al] ( |2007[ ). 



Stellar flux and precision estimator as a function of Wav. Setting 



We found out that the detection and superposition thresholds 
were the factors that decided between the two domains. For val- 
ues bigger than 5%, H band zone dominates. Figure [T] depicts 
the output of this simulation, with the spectral content increas- 
ing along the y axis. Over-plotted is the squared root of the av- 
erage stellar flux, in arbitrary units as well. The final precision 
depends on the product of these two factors and the objective is 
to maximize it. Note the presence of two bumps representing the 
wavelength domains where atmospheric transmission becomes 
so low that the the precision estimators fall along with stellar 
flux. 

The main qualitative difference between the two domains in 
what concerns spectral information can be summarized as: while 
in the K band the spectral information is extracted from more nu- 
merous and shallower atmospheric lines, in the H band it is ex- 
tracted from less numerous, deeper lines. The crowding of lines 
in the K band increases the risk of mutual pollution between the 
two sets of lines (telluric and stellar) which makes precise mea- 
surements more difficult in K band than in H band. 

The final choice was the atmospheric transmission window 
in the H band, and in particular the setting (36/1/i). This setting 
provides for deep, sharp atmospheric lines of CO2 in detectors 
1, 2 and 4. A special attention was given to image the highest 




Fig. 1. Precision as a function of wavelength setting, in arbitrary 
units, as described in Sect. 3.2 {histogram, red in the electronic 
version). Over-plotted is the squared root of the average solar 
flux, as seen by CRIRES in a given wavelength setting, in arbi- 
trary units, (thick line, green in the electronic version). All set- 
ting short-wards of 2.4 fjm are depicted. H band favorable zone is 
located between settings 60-75 and K band's around 110- 115. 
Note that settings where numbered as a function of increasing 
wavelength reference, following the order presented in CRIRES 
manual. 

number of deep stellar lines in detectors 1 and 2, as the blaze 
function provides for a much higher transmission in these two. 
The detector 1 is completely covered by a comb of 22 lines, de- 
tector 2 has 14 lines on it, and detector 4 images 20 lines. Deep 
lines, with relative depths of 10-50%, provide for an homoge- 
neous complete coverage of the first detector; in the second one 
the depths range from 40 to 10% for the first two thirds of the 
detector and are smaller than 10% for the last third; finally, lines 
with a relative depth of 10-55% cover 90% of the fourth detec- 
tor. It is very important to note that there are no other absorbers 
than CO2, providing for both an unambiguous identification of 
spectral features and a well defined continuum around the sharp 
reference lines. For illustration purposes, in Fig.|2]we plot the at- 
mospheric transmission measured by observing the telluric stan- 
dard HD 96 146. 

The simplified study presented here does not intend to have 
the last word on choice setting on CRIRES for RV studies. The 
optimal setting depends not only on the aforementioned super- 
position parameters chosen but also on the stellar spectra, the 
mean RV of the target and even the line identification algorithms 
chosen to pinpoint the lines present in each spectra. More gener- 
ally than that, it depends on the RV determination method used. 
Our objective here is to provide for spectrally separated stellar 
and atmospheric lines, so that the spectral information carried by 
these well-defined lines can be fully exploited. 

As a concluding remark on the issue we would like to note 
that the choice of employing the atmosphere as a wavelength 
reference was prompted by very practical issues. At the time of 
our observations there was no CRIRES gas cell available capable 
of providing deep absorption lines short-wards of 2.5 fim. 

An entrance slit of 0.2" was chosen, allowing for a resolu- 
tion of ~ 100 000. In order to reduce guiding errors (expected 
to have 1-1 translation on the lines profile and thus on the mea- 
sured RV) two cautionary measures were taken. First of all, no 
AO was used, in order to make sure that the star PSF was spread 
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Fig. 2. Atmospheric transmission at the chosen observing set- 
ting, for the detectors 1 (top), 2 (middle), and 4 (bottom). These 
spectra were obtained by co-adding and then normalizing the 
telluric standard exposures taken during the first night. 



and larger than the slit width. On top of that, observations during 
bad seeing periods were requested. The objective was to provide 
enough photons to the slit viewer/guiding camera to ensure a 
good centering of the star on the slit, and reduce photo-center 
effects on RV. While these effects affect both star and reference 
in the same way, minimizing their cause is the best way of en- 
suring a precise correspondence between the two. 



3.3. Observations 

The TW Hya campaign was performed under Director 
Discretionary Time (DDT), with program number 280.C- 
5064(A). The extensively studied TW Hya star was observed five 
times, from the 22 nd to the 28''' of February 2008, promptly fol- 
lowed by a RV velocity standard, HD 108309. This star has been 
observed by HARPS high-precision programs and is known to 
be stable down to 1 . 1 m/s. Between the science target and the RV 
standard the instrumental set-up remained fixed; no instrumen- 
tal changes were allowed for. The objective was to track possible 
zero-point drifts of our wavelength calibration, which could have 
both instrumental or atmospheric origin. The observations were 
made using the wav. setting (36/1/i), characterized by the ref- 
erence wavelength 1592.7 nm (defined as the wavelength at the 
center of detector 3). 

Originally, a set of 10 RV measurements of TW Hya on 
preferably consecutive nights was requested. However, due to 
scheduling difficulties aggravated by multiple instrument inter- 
ventions, the period ended with only 5 out of the 10 measure- 
ments being taken. 

CRIRES calibration plan comprises an RV monitoring of 
radial velocity standard stars with the objective of studying 
CRIRES' stability. From this data we recovered the observations 
on Gl 86. The observations were made in the 5 th , 9 , IT 1 ' and 
18''' of January 2008. The observations were made with a slit of 
0.3", in the setting (36/1/n), and no AO was used. 

Table[T]summarizes the properties of the observations. 

4. Data Reduction & Analysis 

Data were redu ced using a custom pipeline written in IRAF's 
Clf1( Tody|l993] l. The procedures were tailor suited for CRIRES 
spectra using first the science verification data and later the data 
sets analyzed in this paper. It provides for automated dark and 
non-linearity corrections (using the nonlinearity coefficients pro- 
vided by ESO), as well as flagging and replacement of bad pix- 
els. The images are corrected from sensitivity variations by di- 
viding by a flat-fielding corrected from blaze function effect. The 
nodding pairs are mutually subtracted and the orders' tracing 
is fitted using cubic splines. An optimal extraction algorithm is 
employed for extracting the spectra ( |Horne 1 986 ). The objective 
of creating a dedicated pipeline is not only to optimize IRAF 
tunable knobs to provide for an optimized reduction but also to 
ensure a consistent and homogeneous reduction for all spectra. 

The computation of RV from extracted spectra was per- 
formed using a Python and FORTRAN pipeline heavily based 
on HARPS' DRS (Data Reduction Software). The CQ 2 lines' 
central wavelength is provided by HITRAN database (Rothman 
|et al.|200"5j ), with an accuracy of 5 - 50 m/s. The high S/N of the 
spectra and the wide separation between spectral lines allow us 
to fit telluric lines easily and determine their centers. For each 
frame, the wavelength solution was obtained by fitting by a third 
degree polynomial that matches the fitted centers with the theo- 
retical wavelengths. Earth's RV in its orbit around the Sun is cal- 
culated using the ephemerides of Bretagnon & Francou (1988), 
which allows the calculus of barycentric julian dates as well. By 
using the Earth's RV one can transpose the wavelength solution 
to the Solar System barycenter. RVs are then calculated by cross- 
correlating the spectra against a template mask as described in 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Object 


DIT [s] 


# of Nodding pairs 


Reference Wavelength [nm] 


slit width ["] 


Set-up comments 


RVstd 


45 


2 


1592.7 


0.2 


No AO & No Gas-Cell 


TW Hya 


180 


2 


1592.7 


0.2 


No AO & No Gas-Cell 


Gl 86 


30 


3 


1588.4 


0.3 


No AO & No Gas-Cell 



Table 1. The summary of the observation details on the RV std, TW Hya, and Gl 86. All NDIT were set to 1. 



Baranne et al. ( 1996[ l. In our case we use weighted masks (Pepe 
et al.||2002 1 to take into account the impact of photon noise on 



Composite Spectrum on DET1 



lines of different depth. We created two sets of masks: one from 
PHOENIX stellar atmosphere models (Barman et al. 2005} and 
other by measuring the position and depth of each stellar line di- 
rectly on the co-added reduced spectra. The center of the cross- 
correlation function is determined by fitting a Gaussian function. 
We note that we have then the same line fit for both our reference 
and our science target lines. IP variations, which affect the two 
line sets in the same way, are fitted by the same procedure and 
its impact on RVs cancels out. 

As described in Sec. 2, extreme caution must be taken in or- 
der to reduce mutual pollution between stellar and atmospheric 
spectra to a minimum. By knowing earth's RV and the approx- 
imate RV of the targej^] one can calculate the approximate po- 
sition of each telluric and stellar line and discard those that su- 
perpose themselves. While this might lead to a significant loss 
of lines, it greatly reduces the systematic error introduced by 
blended lines. As a result the photon noise calculations on the 
remaining lines provide for a realistic assessment of the spectral 
information carried and achievable precision. 

It is important to note that the RV is computed for each nod- 
ding position spectrum and detector. The RV determination is 
thus independent for each frame, as the always present reference 
lines are fitted independently, and their centers are used in the 
construction of a dedicated wavelength solution. The dispersion 
of the RV within the nodding cycle allows then for the direct cal- 
culation of the uncertainty on the averaged value for the nodding 
cycle. The error bar is simply the r.m.s. of the independent nod- 
ding measurements divided by the square root of the number of 
measurements. 

Along with each RV calculation, the bisector inverse slope 
was computed on the stellar cross-correlation function following 
the procedure described in Queloz et al. ( 2001[ ). 



5. Results 

Due to the very steep blaze function, the low S/N of TW Hya 
spectra collected in detector 4 made atmospheric line fitting im- 
possible. On the other hand, the low number of unblended stellar 
lines in detector 2 combined with the lower count number led to 
a very high photon noise on the CCF function. For the sake of 
consistence and to provide for a simple and robust comparison, 
the standard RV presented were calculated only using the first 
detector as well. In Fig. [3] we plot the spectra of HD 108309, 
on detector 1, for illustration purposes. Out of the 22 lines, 14 
unblended lines were used to calculate the wavelength solution; 
these lines have an average depth of 35. 1 ±12.8% and FWHM = 
3.11±0.31pxl. 

All RV values presented here were obtained with a template 
mask built from the reduced CRIRES spectra. Results derived 



3 Note that the RV variations induced by putative planets are of the 
order of 100 m/s and thus much smaller than the line's FWHM. As a 
consequence, these small amplitude fluctuations can be ignored at this 
stage. 




Fig. 3. Spectrum of RV standard, HD 108309, as seen on detector 
1 . These spectra are the result of the stacking and normalization 
of the 4 nodding spectra obtained during the first night. 



from masks built based on PHOENIX models yielded a lower 
precision. The dispersion of the RV standard measurements, the 
most reliable precision indicator among our measurements, is of 
~25 m/s for PHOENIX masks. The reason for this difference is 
not known. 

The RV measurements of TW Hya and the associated stan- 
dard star are presented in Tab.|2]and plotted in Fig. [4] The disper- 
sion of the standards' RV measurements is compatible with the 
error bars, while for TW Hya that is not the case. Simultaneous 
optical RV determination were performed by CORALIE to back 
CRIRES observations, already described in detail in Huelamo 



et al.| ( |2008] >. TW Hya IR RVs are not compatible with the or- 
bit drawn from optical RV measurements, with an amplitude of 
roughly 200 m/s, defined using CORALIE data and reproducing 
perfectly Setiawan et al.l (|2008[l data. To test the hypothesis of 



a smaller RV semi-amplitude in the IR we let this quantity and 
the systemic velocity as free parameters, fixing all others. The fit 
delivers a semi-amplitude of 80.50±6.83 m/s with a r.m.s. (O-C) 
around the fit of 7.93 m/s. The fit is over-plotted on the data. 

The RV measurements of Gl 86 are presented in Tab. [3] and 
Fig [5] The orbit presented in Queloz et al. ( 2000 ) is superposed 
on the data. At first, only the systemic velocity and the To were 
allowed to vary. The reason for letting To as a free parameters 
comes from the fact that the error propagation from the 2000 or- 
bit leads to a very large error on this parameter, as of 2008. The 
uncertainty adds to almost 8 days, which is roughly half of the 
orbital period. The agreement between the published curve and 
CRIRES measurements is shown by the small dispersion around 
it, of only 5.4 m/s. Alternatively, one can fix To and let P to vary; 
the P remains unchanged within error bars, and the dispersion 
of the residuals remains roughly the same. We recently obtained 
several CORALIE (Baranne et al. 1996) measurements on this 
star which, together with the 2000 data, allowed us to determine 
the To in an independent way. This procedure delivered a value 
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RVstd MJD [d] RV[km/s] Ph. Noise [m/s] TWHya MJD [d] RV [km/s] Ph. Noise [m/s] 



2454519.805826 


30.3699 


10.5 


2454519.795056 


12.1094 


20.6 


2454519.806647 


30.3600 


10.1 


2454519.797602 


12.0837 


19.1 


2454519.807307 


30.3452 


10.1 


2454519.799998 


12.1662 


19.4 


2454519.808129 


30.2390 


10.7 


2454519.802579 


12.0965 


18.9 


2454522.666534 


30.3146 


16.8 


2454522.655881 


12.2158 


32.9 


2454522.667368 


30.3778 


16.4 


2454522.658416 


12.1857 


26.0 


2454522.668039 


30.3330 


16.2 


2454522.660823 


12.2122 


29.3 


2454522.668861 


30.3402 


16.7 


2454522.663370 


12.2435 


29.9 


2454523.728265 


30.2983 


14.4 


2454523.716034 


12.0599 


26.1 


2454523.729098 


30.3300 


14.7 


2454523.718569 


12.0379 


23.8 


2454523.729770 


30.3419 


16.4 


2454523.720965 


12.1981 


25.9 


2454523.730592 


30.3281 


15.7 


2454523.723534 


12.0555 


25.8 


2454524.749738 


30.3288 


13.4 


2454524.734402 


12.0436 


24.6 


2454524.750571 


30.3301 


12.3 


2454524.736948 


12.1082 


22.8 


2454524.751242 


30.3526 


11.6 


2454524.739344 


12.0788 


34.0 


2454524.752064 


30.3389 


12.5 


2454524.742249 


12.0935 


26.1 


2454525.733788 


30.3323 


11.5 


2454525.713207 


12.2017 


22.2 


2454525.734610 


30.3690 


10.0 


2454525.715765 


12.1959 


18.0 


2454525.735270 


30.3403 


9.6 


2454525.718161 


12.1969 


19.3 


2454525.736092 


30.3208 


9.8 


2454525.720730 


12.1501 


19.7 



Table 2. The RV calculated per nodding position for the RV standard, HD 108309 and TW Hya, and estimated photon noise of each 
exposure. 



TWHya 




0.0 1.0 2.0 3.0 4.0 5.0 6.0 

JD - 2454520.0 [days] 



Fig. 4. RV variation of the standard star {left) and TW Hya {right) as a function of time. The error bars are drawn from the dispersion 
of the N=4 nodding positions, divided by 2 ( VA0, as described in Section 4. The optical orbit, with a different K (as described in 
Section 5) is over-plotted on the points. 



compatible, within error bars, with the one provided by CRIRES 
points alone. It is important to note that there is a measurable 
drift due to the presence of an outer companion (Eggenberger 
|et al.|[2003] l and the determination of the To at the moment of 
CRIRES data (when there are no CORALIE points) is affected 
by the shape of the unknown orbit. In any case, this To value 
can be used to fit an orbit in the CRIRES data, leaving only the 
systemic velocity as the free parameter. This delivers a disper- 
sion of 16 m/s, which, while being bigger, shows the agreement 
between the orbit and CRIRES data. 

We can estimate the precision expected on computing the 
cross-correlation function on each spectra, drawn from photon 
noise alone. To do so we use the method presented in Bouchy 
|et al.| ( |20Q 1 ] >. We considered the particular case of a gaussian 
absorption line and applied it to the stellar lines present in the 
correlation mask (the unblended ones actually used in the cor- 



relation, and of interest to us). This yields the precision on each 
individual spectra, i.e., on each nodding position. To calculate 
the error on the averaged night value we need to divide it by the 
square root of the number of independent nodding positions. 

We can then compare three different estimators of the 
reached precision: 

- the weighted dispersion of the night-averaged RV points - 
External Dispersion; 

- the error on the averaged night value, given by dispersion on 
the nodding cycle RV values divided by V^V - Intra-Night 
Dispersion; 

- the average photon error estimator value (as described 
above) divided by V^V - Photon Noise. 

Table|4]presents all these results. Note that these errors were 
estimated in such a way that they can be easily compared, as they 
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External Dispersion [m/s] 


Intra-Night Dispersion [m/s] 


Photon Noise [m/s] 


(O-C) [m/s] 


RVstd 


5.77 


7.03 


6.48 




TWHya 


54.57 


12.12 


12.10 


7.93 


Gl 86 


122.47 


12.77 


7.62 


5.41 



Table 4. The different RV precision indicators on the RV std, TW Hya, and Gl 86. The dispersion around the described fits are 
presented for TW Hya and Gl 86. 



Gl 86 MJD [d] 


RV [km/s] 


Ph. Noise [m/s] 


24544 /0.5 190 lb 


55.0081 


1 1 
18.1 


24544/0.519/13 


55.595 / 


It). / 


24544 /0.520 1 99 


55.6/05 


1 O 1 

18.3 


H A £ A AH(\ ZHt\Q AH 

24544 /0. 52084 / 


55.661 / 


on a 


H A £ A AH(\ ?T1 1AQ 

24544 /0.52 1 345 


55.0/03 


1 1.0 


H A £ A AH(\ CI 1 fl£n 

24544 /0.52 1 969 


rr Tin 

55. 1 ill 


1 n 
1 /.8 


HAZAAHA ZAH^iAH 

24544 /4.54 /24 / 


56.3486 


"1 
8. / 


HAZAAHA C/iTmn 

24544 /4.54 /930 


56.2390 


Q A 

8.4 


2454474.548508 


56.2105 


8.7 


2454474.549145 


56.1806 


8.5 


2454474.551020 


56.1226 


8.5 


2454474.551656 


56.3237 


8.1 


2454482.535667 


55.7555 


12.9 


2454482.536315 


55.7780 


11.6 


2454482.536813 


55.7990 


11.9 


2454482.537461 


55.7627 


14.3 


2454482.537970 


55.7263 


13.6 


2454482.538618 


55.7589 


12.0 


2454483.517437 


55.8188 


11.9 


2454483.518085 


55.5305 


11.8 


2454483.518571 


55.6780 


12.2 


2454483.519207 


55.7346 


12.2 


2454483.519694 


55.4173 


13.3 


2454483.520342 


55.6021 


11.5 



cr [m/s] P [%] 



Table 3. The RV per nodding position for Gl 86, and estimated 
photon noise associated. 




-0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 



Fig. 5. Phase-Folded plot of Gl 86 Rv measurements. The error 
bars presented are drawn from the dispersion of the six nodding 
positions, divided by V6, as discussed in Section 4. 



provide for alternative ways of looking at the same quantity. In 
the cases of TW Hya and Gl 86, these can be compared with the 
dispersion around the fit. 



5 
6 
7 
8 
9 
10 



84.1 
66.3 
49.9 
36.6 
26.6 
19.5 



Table 5. The probability of obtaining a dispersion of 5.77 m/s 
on a set of 5 points characterized by dispersions in the range 
5-10 m/s. 



6. Discussion 

In our previous work, ( Huelamo et al.|2008]|, we reached a pre- 



cision of 35 m/s over a time span of roughly 1 week. This value 
was obtained in relative RV: science target - RV standard. This 
procedure allowed for a correction of the systematics still present 
in the data after our first reduction. Seifahrt & Kaufi (2008), a 
contemporary study, reached a precision of 20 m/s on one night. 
A dedicated effort was put in order to track systematics and un- 
derstand their possible causes. As a result we could now improve 
the analysis procedure in order to extract the RV in an optimized 
way, leading to a significant increase in precision. We have no 
need to calculate relative RV in order to correct for systemat- 
ics anymore. The RV standard dispersion is now close to photon 
noise, showing that there are no systematics present down to 5- 
10 m/s. The same is valid for TW Hya and Gl 86 (O-C) residuals 
which show a similar dispersion. 

The RV standard displays a variation of less than 6 m/s 
(r.m.s.) over 1 week. The probability that a dispersion of 5.77 m/s 
or lower is obtained on 5 points belonging to a distribution 
with <x=6.48m/s (the estimated photon noise) is of 58.2%. If 
we consider 5 points drawn from a distribution characterized by 
cr=7.66 m/s (the measured internal dispersion), the probability 
drops then to 40.8%, but it is still high. We calculated the proba- 
bility of drawing a dispersion of 5.77 m/s from distributions with 
cr from 5 to 10 m/s. The results, presented in Table [5] show that 
we cannot exclude that the intrinsic dispersion of our distribution 
is slightly larger than the measured one. As a result we consider 
the attained precision, in a conservative way, as being of 6-8 m/s. 

It is interesting to compare the values of intra-night disper- 
sion with photon noise estimations, for the three targets. When 
measuring the intra-night dispersion, one takes into account not 
only the photon noise but also the effect of systematic errors on 
RV. A straightforward example of such are the errors induced 
by nodding and locating the target at different positions in the 
slit, known to have two slit jaws that are not perfectly parallel. 
While the values are very similar for the RV standard and TW 
Hya, that is not the case for Gl 86. A possible cause is the usage 
of a larger entrance slit, and the consequent introduction of ad- 
ditional sources of error, such as stellar photo-center variations. 

The RV values for the standard and TWHya are compat- 
ible, within error bars, with the previously published values. 
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This is mostly due to the comparatively large dispersion an- 
nounced previously (35 m/s) and error bars with even larger am- 
plitude (30-60 m/s). However, now, given the higher precision 
of the measurements, a low-amplitude signal emerges from the 
noise. As previously discussed in (Huelamo et al. 20081, a lower- 



amplitude signal in the IR than in the optical is expected if a spot 
is responsible for the RV variations instead of a planet. The pre- 
viously presented tests stand as well. One can calculate the orbit 
parameters using optical observations. The probability that the 
IR RV points are drawn from this orbit by adding an error com- 
parable to the instrumental precision is below lxl0~ 3 %, show- 
ing that the spot theory prevails. 

The RV determinations for Gl 86 match the orbit as estab- 
lished by the announcement paper. It is important to stress that 
for these fit the semi-amplitude, P, e and a) were fixed on the 
published values in |Queloz et al.| ( |2000| l. This shows that our 
method yields the correct value for RV signals, reproducing the 
bona-fi.de planet detected at optical wavelengths. 

Given the very small number of points for both TW Hya and 
Gl 86, the dispersion around the respective orbits is to be taken 
with caution. At this level the number of free parameters ap- 
proaches the number of points and the noise contribution is fitted 
as well, leading to a^ 2 < 1. As a representation of this, one can 
use the TWHya data in a slightly different way. One may use 
the RV values obtained for each nodding position, and assign as 
error bars for each given by the photon noise estimation. The fit 
delivers, within error bars, the same orbit parameters (K and y), 
a dispersion of the residuals of 35 m/s and a;^ 2 =1.63, showing 
that the error bars are underestimated by photon noise alone. At 
this point our objective is merely to show that the RV points are 
drawn from the orbits. They do not define the orbits and the dis- 
persion around the same is not representative of the true errors 
for such a low number of points. 

A question that arises naturally is "What are the factors lim- 
iting the RV precision delivered by this method?". This is diffi- 
cult to evaluate mostly because the data set we have in hands is 
small. However, on theoretical ground, some educated guesses 
can be made. The precision on each of the telluric lines central 
wavelength is not quantified on HITRAN database, from which 
the values were drawn. The lines are merely flagged as having 
a precision between 5 and 50 m/s. It is thus not surprising that 
a precision of 5 m/s can be reached in a whole, but this might 
become the limiting factor if one pushes the S/N to much higher 
values. The presence of a systematic wavelength assignment er- 
ror can also be a trivial source of error. On the other hand, at- 
mospheric phenomena such as localized jet streams may imprint 
an RV on the reference lines themselves. However, its is hard to 
evaluate how large these effects are when the we integrate along 
the line of sight, up to the outer atmosphere. The conjugation of 
three factors point to CO2 features as some of the most stable of 
the atmosphere, though: 1) the abundance of CO2, 2) the con- 
stance of its volume mixing ratio of ~300 ppmv up to 80 km, 3) 
the low absorption line intensities. When taken together, these 
properties reduce the impact of localized air mass movements 
on the RV. With the current set of data we can only evaluate the 
short-term stability of the atmospheric lines. An extended data 
set with observations under very different atmospheric condi- 
tions is now under evaluation. The results will be the subject of 
a forthcoming paper (Figueira et al. in prep.). 

Temperature and pressure changes are known to be one of 
the main causes for spectrograph's IP variations. It is impor- 
tant to note that CRIRES is cryogenic and thus stabilized in 
temperature. At the time of writing, the thermal stability of the 
most important components of the spectrograph is of ~0.05 K for 



the prism, 0.2 K for the grating and 0.3 K for the Three-Mirror- 
Astigmat. As a consequence, the IP remains constant and its im- 
pact on RV measurements, by introducing line asymmetries, is 
reduced. Another possible lines profile variation comes from at- 
mospheric phenomena. In what concerns these, and as referred 
before, the eventual variations are seen both by our reference 
and target spectra. The drawback is that this shows that typi- 
cal bisector analysis cannot be used in a straightforward way in 
stellar spectra, as one must remove the effect of the variation of 
atmospheric lines' bisector before. It is then not a surprise that 
bisector shows no correlation with the RV. Still, the low num- 
ber of points and the possibility of atmospheric profile variations 
makes this diagnosis very weak. The data do not provide for a 
finer analysis, which is not, anyway, within the reach of this pa- 
per. 

The contribution of photocenter errors was negligible in our 
error budget but that might be not the case if the seeing becomes 



very low, as experienced by Bouchy et al. (2005). However, our 
situation is fundamentally different, because our reference lines 
will be affected by affected by these photocenter movements in 
the same way than the object lines. In any case, a simple way of 
reducing these effect is the usage of a fiber entrance, providing 
for an homogeneous illumination over a wide range of seeing 
conditions. 

The fact that PHOENIX-derived masks led to a lower pre- 
cision than masks derived from the data themselves is quite in- 
sightful on its own. It states that the line's position precision as 
calculated by state of-the art models is lower than the one that 
can be measured from high S/N spectra. Some RV determination 
methods are more sensitive to this fact than others. In particular 
those that obtain the RV by searching for a match between stellar 
models and observed spectra struggle against an additional and 
non-negligible source of noise. 

While the application of gas cells on general-purpose slit 
spectrographs such as CRIRES is very tempting, their yielded 
precision remains to be assessed. As already discussed by 
Mahadevan & Ge ( 2009[ ), the usage of the gas cell as a simul- 
taneous reference superimposed on the stellar spectrum is lim- 
ited by the definition of the stellar continuum. In the case of M 
stars, one of the main drivers for RV in the IR, the continuum is 
shaped by a dense forest of unresolved spectral lines. These fea- 
tures are convolved with the atmospheric transmission, which is 
variable through time. These factors together make RV determi- 
nations much more demanding, time and photon consuming and, 
in many cases, degenerated. In our particular approach, such a 
poor continuum definition would introduce systematic errors in 
the absorption lines wavelength assignment, which would prop- 
agate to the final calculated RV. 

Several infrared spectrographs dedicated to RV measure- 
ments are being designed, such as NAHUAL (Martin et al. 2005 ) 
or SPIROUQ These spectrographs aim at covering a wide spec- 
tral range simultaneously, from around 1 to 2.2 pm. If a gas cell 
on the main optical path is chosen as the calibration procedure, 
the yielded precision will vary locally as a function of stellar 
continuum definition and atmospheric subtraction quality. On 
the other hand, using this gas-cell on a twin optical path with 
an homogeneous continuum source or other methods such as a 
laser coupled to a Fabry-Perot resonant cavity, allows a simple 
workaround of these problems. However, as discussed before, 
this approach is only possible in a fiber-fed, stable spectrograph. 
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7. Conclusions 

We proved that RV measurements with a precision of 5-10 m/s 
are within reach of CRIRES while using carefully selected atmo- 
spheric features as a wavelength reference. This precision was 
obtained on an RV standard, HD 108309 over a time span of 
roughly one week. TW Hya RV variation in the IR is not compat- 
ible with its optical counterpart, showing that a stellar spot is the 
best explanation for the observed RV variations. This confirms 
our previously published results (Huelamo et al. 2008) with bet- 
ter precision and consistency. We also showed that the RV vari- 
ation of a star harboring a planet, Gl 86, is well reproduced by 
CRIRES data with a precision which is compatible with photon- 
noise estimations and at the same noise level. 

The usage of the Earth's atmosphere as a free gas-cell is still 
to be fully assessed, and the outcome of these studies will, in 
all probability, shape the design of future dedicated IR spec- 
trographs. In order to do so more lengthy data sets spanning 
a longer time-span are required. However, while this method 
seems promising and able to compete with regular gas-cells ap- 
proaches, it remains to be proven that a higher precision is pos- 
sible. Strong limitations might arise below 5 m/s steaming from 
sources such as the knowledge of the precise wavelength of tel- 
luric lines or their stability over long time-scales. 
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